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Abstract 

We study the dipole anisotropy in the arrival directions of high energy CR elec- 
trons and positrons (CRE) of Dark Matter (DM) origin. We show that this 
quantity is very weakly model dependent and offers a viable criterion to dis- 
criminate among CRE from DM or from local discrete sources, like e.g. pulsars. 
In particular, we find that the maximum anisotropy which DM can provide is 
to a very good approximation a universal quantity and, as a consequence, if a 
larger anisotropy is detected, this would constitute a strong evidence for the 
presence of astrophysical local discrete CRE sources, whose anisotropy, instead, 
can be naturally larger than the DM upper limit. We further find that the 
main source of anisotropy from DM is given by the fluctuation in the number 
density of DM sub-structures in the vicinity of the observer and we thus devote 
special attention to the study of the variance in the sub-structures realization 
implementing a dedicated Montecarlo simulation. Such scenarios will be probed 
in the next years by Fermi-LAT, providing new hints, or constraints, about the 
nature of DM. 



1. Introduction 

High energy Cosmic Ray (CR) positrons are promising targets for indirect 
searches of Galactic particle Dark Matter (DM) [T] . The recent results reported 
by the PAMELA 2J and Fermi collaborations [31 |U [S] on the positron fraction 
e+/(e+ + e^) and on the (e+ + er) (CRE) spectra in the GeV TcV energy 
range show large discrepancies with standard astrophysical model predictions 
and have indeed raised a large number of interpretations in terms of DM. In 
particular, it has been shown that a good fit of the PAMEL A/Fermi data can be 
achieved with e^e~ produced by a DM particle of '^TeV mass annihilating or 
decaying predominantly via leptonic channels, and several models realizing this 
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scenario have been proposed [HI HI UHl E] • However, interpretations based on 
discrete astrophysical extra sources (like e.g. pulsars, or stochastic local sources) 
[lIl[12l[T3lllllll5l[16l[17l[18] have been shown to provide equally good fits to 
the data (see e.g. [IH] for an extended critical review of the subject). 

It is, however, very unlikely to distinguish the two scenarios using as observ- 
ables only the CRE fluxes, even with the larger statistics expected in the future 
[20] . It is thus mandatory to find other observables accessible to experiments, 
that are as much model independent as possible and can provide a clear discrim- 
ination between a DM dominated scenario and an "astrophysically" dominated 
one. The intrinsic degree of dipole anisotropy in the arrival directions of high 
energy CREs expected from a DM scenario, 5dm, is indeed insensitive to many 
uncertainties, and constitutes, to a good approximation, a universal character- 
istics of galactic DM, as we will show in the following. The reason why the 
dipole anisotropy has a very weak dependence on the many unknowns involved 
in the problem is, on the one hand, the very short 1 kpc) electron path above 
~ 100 GeV which makes this quantity very local in origin, on the other hand, 
the fact that it is a flux ratio (see Eq. [T]) so that most of the uncertainties cancel 
each other in the ratio. 

Furthermore, the anisotropy signal from DM is intrinsically very different 
from the one due to local discrete sources. Pulsars are rare (their number in 
the Galaxy is estimated to be around 10^ [H]) and powerful objects and can 
induce very large anisotropics typically dominated by a single or a few nearby 
objects. On the other hand, the number of galactic DM substructures is 0(10"'^'') 
and they produce a "collective" anisotropy which is never dominated by a single 
clump. The flux from a very nearby clump would be always accompanied by the 
large, dominant and almost isotropic flux from the whole population of clumps, 
which washes out the single clump anisotropy. Therefore, the dipole anisotropy 
offers a viable criterion to discriminate among CREs produced by DM or in 
local sources. 

Anisotropies in the DM component have been also studied in gamma-rays 
(see for example [HI [531 (24 ) which have clearly the advantage of being inde- 
pendent of the choice of a diffusion models. On the other hand, as we will show, 
CREs anisotropies are also fairly independent of the propagation setup over a 
wide range of possible diffusion models and thus offer an interesting comple- 
mentary anisotropy probe of DM. 

Besides the anisotropy from DM and local discrete astrophysical sources, 
there is also a third source of anisotropy which needs to be considered in order 
to have a complete picture, i.e. the anisotropy from the large scale distribution 
of the astrophysical sources considered as a whole. This component, as we 
will show, gives generally a smaller anisotropy with respect to the first two 
components above. 

On the experimental side, Fermi-LAT recently placed the first upper limits 
on the integrated dipole anisotropy of the arrival directions of CRE with E > 
60 GeV [55] , and there are prospects for its actual observation after a few years 
of data taking, if local pulsars contribute significantly to the CRE fluxes above 

100 GeV [H]. Also AMS-02 [26] is now taking data, but its sensitivity 
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to CRE anisotropy is much lower than the one of Fermi [50]. We will show 
that if Fermi- LAT or future experiments will find an anisotropy larger than the 
maximum DM anisotropy we derive here, then a dominant DM contribution 
to the CRE anisotropy can be excluded in a basically model independent way, 
pointing instead to local discrete astrophysical CRE sources as the main source 
of anisotropy. Therefore, the observation of anisotropy at the level within reach 
by Fermi-LAT in the next years will be able to constrain significantly the flux 
of CREs possibly contributed by DM annihilations in the Galaxy. On the other 
hand, to identify pulsars as responsible for a possible anisotropy would require 
at least a careful analysis of their spectral characteristics and of the direction 
and intensity of the anisotropy. 

This paper is structured as follows: in Section [2] we describe how we com- 
pute the interstellar CRE density due to astrophysical sources and due to DM 
annihilations in the smooth halo and in substructures. In section [3] we detail 
how we simulate the distribution of galactic DM substructures. In section |4] we 
compute the total intrinsic DM anisotropy (i.e. the anisotropy when, ideally, 
only DM contributes to the total CRE emission) while in sectionjsjwe discuss the 
anisotropy for a mixed scenario in which both DM and standard (non discrete) 
astrophysical sources contribute to the CRE flux. Section [6] is finally devoted 
to our final comments and conclusions. 

2. DM intrinsic electron anisotropy 

In the diffusive approach, the dipole anisotropy can be written as |27j 

where D is the diffusion coefficient, /3c and are the CRE velocity and flux 
respectively. The total DM contribution to the e+e~ fluxes can in general be 
written as the sum of two components, (pDM = (f'h + (f's, where 4>h is the contri- 
bution from the smooth halo while 0^ is the contribution from the substructures. 
For each contribution, we have 

where G is the Green function associated to the transport equation [27], is 
the DM density at the Solar System position and dN^/dE' is the annihilation 
spectrum into e+e^. The term p''/^ [x'Y is defined as {ph{x') / pof in the case 
of the DM halo density [i — h), while in the case of the substructures {i = 
s) is written as pl^^ [x'Y = iPji^')/ Pq) : with the sum running over the 
substructures and pj representing the DM density of the single substructure. 

The large scale anisotropy from astrophysical sources will be estimated from 
the interstellar electron density computed numerically with the DRAGON code 
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pS] . We will indicate in the following this contribution with the acronym AP. 
It is worth stressing, however, that this does not include the anisotropy from 
discrete sources like pulsars which is in general larger (see [15] and [16]) and 
which is the contribution that we want to isolate in the case of detection of an 
anisotropy larger than the one provided by DM. 

2.1. Electron propagation 

We solve in the stationary limit dN/dt — the well known diffusion-loss 
equation [37] 

^ - V • (d{E)Vn) - A AT) = Q^E, x) , (3) 

where N is the particle number density, b{E) represents energy losses, D{E) = 
Dq{E/3 GeV)" is the (spatially constant) isotropic diffusion coefficient and 
Q is the source term. Since \5\ is measured at i? > 60 GeV only diffu- 
sion and continuous energy losses affect significantly the propagated spectra, 
so we can neglect reacceleration and convection. Moreover, at high energies 
leptons cannot travel more than a few kpc [27], hence we assume b{E) = 
1.6 X 10~^^(i?/GeV)^ GeV/s corresponding to a magnetic field and interstel- 
lar radiation field constant over the relevant propagation region, whose vertical 
height scale we fix as i = 4 kpc. As a further consequence, the effect of bound- 
ary conditions at > 60 GeV is negligible. It can be checked however that 
changing L in the range allowed by CR nuclei constraints [29] does not produce 
a significant effect. Given that the largest effect on the dipole anisotropy is 
expected to come from the rigidity dependence of the diffusion coefficient, we 
consider three different models of diffusion: Kolmogorov-like turbulence (KOL) 
with a — 0.33 and Dq = 5.8 x 10^* cm^s~^, Kraichnan-hke turbulence (KRA) 
with a = 0.5 and Dq = 3 x 10^^ cm^s"-^, and a last one (HA) with a = 0.7 and 
Do = 1.28 X 10^8 cm^s"^ The values of Dq are in agreement with CR nuclei 
observations 15J. Further details on how we solve the diffusion equation and 
derive the related Green functions are given in [Appendix A[ 

2.2. Dark Matter distribution 

Beside the diffusion setup another ingredient required to derive the DM 
anisotropy is its distribution in the Galaxy. Current highest resolution N-body 
simulations ({30] [31]) find that the DM mass is distributed roughly equally into 
a smooth Halo component and into a further clustered part (DM clumps or 
substructures). With respect to the anisotropy, the halo component gives only 
a mild dipole anisotropy, which for symmetry reasons points toward the Galactic 
Center. The clustered component, instead, gives the main contribution to the 
anisotropy, but is also more difficult to model. We give below the main points 
while full details are reported in [Appendix B[ 

TV-body simulations roughly agree on the mass distribution of substructures, 
predicting a number density scaling like (Via Lactea II [3D|) or rn~j^'^ 
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Figure 1: DM anisotropy as function of energy contributed by substructures in different mass 
ranges for 100 different realizations. The plots refers to the case of a 1 TeV DM particle 
annihilating into fi pairs, for a NFW distribution and a KRA propagation setup. 



(Aquarius ,31]). How substructures are distributed in the smooth halo is how- 
ever more uncertain. We considered the two extreme cases of an unbiased 
distribution where substructures follow the main halo and an anti-biased case 
as suggested by the Via Lactea II simulation [32] . The internal concentration of 
substructures and the effects of tidal disruption are parameterized as in [33] . We 
considered also a very different set of hypotheses (concentration parametrization 
taken from [M] and no tidal effects) finding almost unchanged results, which 
suggests that internal concentration and tidal forces play a minor role on S. 
Finally, we chose a clump mass range 10^^ ~ 10^" Mq, with the upper limit 
coming from constraints due to disk stability [35 . The lower limit is set instead 
following the most common choice in the literature. Our results however do not 
depend critically on the lower limit. For the spatial distribution of the smooth 
component and for the DM distribution inside the substructures we consider 
Navarro-Frenk- White (NFW) [3S] and Burkert [37] profiles. 



3. DM Anisotropy and Clumps Simulations 

With the explicit solution of the diffusion equation and an analytic form 
for the distribution of the clumps it is possible to express the average dipole 
anisotropy resulting from the sum of all the clumps in a completely analytic 
form which we report in [Appendix C[ The result shown in [Appendix C[ is valid 
when the number of clumps contributing to the anisotropy is very large, i.e. for 
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Figure 2: Angle between the direction of the DM dipole anisotropy and that of the Galactic 
Center as function of energy, taking into account substructures in different mass ranges. The 
lowest curve include the contributions from the smooth halo and the entire substructure 
distribution, 10"® IO^OMq, averaged over 100 different realizations. The plot refers to the 
case of a 1 TeV DM particle annihilating into pairs, for a NFW distribution and a KRA 
propagation setup. 



clumps in the small (mc/ <10^ -^0) mass range, because in this case fluctuations 
are much smaller than the average. Clumps of higher mass are instead less 
abundant and the stochastic fluctuations due to the particular realization of 
their spatial distribution can give a sizable contribution to the anisotropy which 
must be taken into account. To this purpose we thus perform explicit Montecarlo 
simulations of the substructures. 

At present, it is computationally prohibitive to simulate the whole pop- 
ulation of the 0(10^^) substructures lying within the diffusive region. We 
thus compute analytically the average contribution from substructures with 
10~® <mci/MQ < 10^, while we compute explicitly the contribution of each 
clump with mci> 10^ Mq. We sampled the distribution of md >10^ Mq sub- 
structures via a MonteCarlo procedure as described in j3H]. We produced 100 
realizations of substructures with 10"'<TOc;/A^o<10^° and 10 realizations with 
lO^<TOci/M0<lQ^. We then averaged our results over the 100 realizations, hav- 
ing checked that the fluctuations induced by the 10 samples with low masses 
are much smaller and can be neglected. In total, we computed the contribution 
of 0(10^°) substructures for each model. 

According to the analytic computation shown in |Appendix C[ the anisotropy 
associated to clumps within a mass decade (10" < itIci/Mq < 10"+^, with 
—6 < n < 2) is independent of the mass decade (i.e. it is independent of n) 
to a very good approximation. Therefore, the anisotropy contributed by small 
clumps results in the sum of several equal contributions. Fluctuations in the 
distribution of the high mass substructures however break the mass scaling and 
lead to large variability of \d\. This effect is more relevant at higher energies, 
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because higher energy electrons probe smaller volumes than lower energy ones, 
thereby being more sensitive to fluctuations in the clump distribution. This is 
confirmed by the results in Fig. [T] (see also |33] for an analogous discussion about 
fluctuations in fluxes) where we show the anisotropy resulting from each of the 
100 simulations. In particular, we plot the anisotropy as a function of energy for 
a DM candidate x with = 1 TeV, annihilating into , assuming NFW 

profile and KRA propagation model, for different clumps mass decades. It is 
clear from the plot that the main effect of the small mass substructures is to set 
the average value of the anisotropy, while high mass clumps {rrici ^ lO'^ -^0) 
are responsible for fluctuations with respect to the average as large as even one 
order of magnitude. We remark that, even if the 100 realizations include only 
clumps in the mass range 10^ < tUci/Mq < 10^°, there are fluctuations also 
for masses 10~^ < uIci/Mq < 10^, which come from the denominator in the 
definition of the dipole anisotropy (that includes the total flux coming from the 
whole clump sample). These fluctuations are small, in agreement with findings 
from [55] . 

Finally, we show in Fig. [2] the direction with respect to the GC of the DM 
dipole anisotropy, averaged over our 100 realizations and for different ranges 
of substructure masses. The plot refers to the case of a 1 TeV DM particle 
annihilating into /i pairs, for a NFW distribution and a KRA propagation setup. 
While the dipole resulting from high mass clumps only can be pointing in a 
random direction, the contribution of smaller substructures stabilizes the dipole 
towards the direction of the Galactic Center. Notice, however, that while at 50 
GeV the dipole points within a few degrees to the GC, at 500 GeV the dipole 
can be pointing up to 25 degrees off. 



4. Intrinsic anisotropy upper limits 

The anisotropy signal that can be measured by our observatories must clearly 
be computed as a suitable composition of the contributions of all the compo- 
nents of CRE fluxes. While we discuss a two-component framework in detail 
in Section [5] we discuss here the single component scenario in which only DM 
contributes to the CRE fluxes above 60 GeV. While this assumption is not 
justified below ~ 200 GeV, where data on the positron fraction [2 [5] imply 
that a possible DM contribution to these fiuxes can be at most of the order of 
30%, the single component scenario becomes more realistic at larger energies, 
where uncertainties on the primary astrophysical CRE component are larger. 
However, we stick to the use of the single component scenario at all energies in 
order to retain the simplicity of the discussion and to emphasize the intrinsic 
properties of the DM anisotropy. For limits below 200 GeV, where a substantial 
astrophysical component is required, this is, nonetheless, a (overly-)conservative 
choice, as discussed in more details in the next section. 

The total anisotropy from the smooth Halo and the substructures can be 
expressed as: 

_ 3D{E) WME) + WME) ... 
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Figure 3: Intrinsic integrated dipole anisotropy of DM, compared to the results for the stan- 
dard astrophysical background in the associated propagation scenario (Sap, black dashed 
Une). Current Fermi upper limits at 95% CL and the sensitivity expected after 10 years of 
data taking (actual limits rescaled by a factor \/T0, accounting for the typical \/N scaling of 
sensitivity with statistics — this represents an optimistic case in which systematics effects do 
not dominate the measurement) are shown only for reference. The points correspond to dif- 
ferent annihilation channels and masses of the DM particle. Both average values and maxima 
are shown. The thick black line shows the upper envelope of the maxima. 
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We consider annihilation in /i, t and quark pairs, for values of the DM mass: 
100, 316, 1000, and 3162 GeV. Q This basically spans all the possible spectra 
which can arise from different particle physics model, apart perhaps the case of 
models where the annihilation is mediated by a light boson as in [TD]. Even in 
this case, however, the spectra are only mildly different from the cases considered 
above and, as we will see, the results are anyway fairly independent of the 
annihilation channel. Figure [3] shows the results for the mean and maximum 
5dm = \Sdm\ for the considered annihilation channels and DM masses, for our 
chosen DM density profiles and propagation setups in the unbiased case. The 
anti-biased case (not shown) gives results smaller by a factor of ~ 5 at 500 
GeV. The maximum over our 100 realizations roughly corresponds to a 99% 
CL. The figures also show the prediction for the anisotropy of the Astrophysical 
Background Sap = \Sap\ calculated with DRAGON for the various propagation 
scenarios considered. As already noted, this does not include the contribution 
to anisotropy from local discrete sources. Figure |4] shows, instead, only the 
envelope of the maxima for the various propagation setups. Both the mean 
and maximum anisotropy increase with energy, as expected from the fact that 
at higher energies smaller and smaller propagation volumes are probed by the 
CREs and the role of fluctuations is more relevant. 

The main result emerging from the above plots is that Sum is almost in- 
dependent of the detailed characteristics of the DM models and distributions 
in substructures and, in this sense, 5dm is a general property of DM. Being 
a ratio, it is very little sensitive to integrated quantities, like the annihilation 
spectrum. Moreover, because CREs propagate only a few kpc distance in the 
Galaxy, Sdm is also little sensitive to the DM spatial profile, in particular on 
whether it is peaked or cored. The anisotropy is also not strongly sensitive to 
the internal concentration of the subhaloes, because diffusion over kpc scales 
smooths out the effect of a possible cusped over-density region. For these rea- 
sons, as we checked, the case of decaying DM gives similar results as the case of 
annihilating DM. Finally, remarkably, also the fluctuations in Sdm are basically 
model independent. Figure |4j in particular, shows that the maximal intrin- 
sic DM anisotropics are almost independent of the diffusion setup, even in the 
extreme case of HA diffusion, and only slightly change for different DM profiles. 

5. Anisotropy upper limits in mixed astrophysical-Dark Matter sce- 
nario 

The results described in Section |4] refer to anisotropy from DM only. We 
now consider again the role of the background. To fix ideas, we will consider for 
the AP background an injection spectrum of CRE dN/dE cx E^^-^^, which we 
then propagate with DRAGON using the same propagation setup as for the DM 
contribution and we normalize such that it accounts for 90% of the CRE flux 



We compute the injection spectra in these models with the numerical package DarkSUSY 
[39] Ho]. 
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Figure 4: Upper limits (UL) on the intrinsic DM anisotropy for all the combinations of DM 
density profiles and propagation models under scrutiny. Also shown for reference are current 
Fermi UL and its expected sensitivity after 10 years data taking. 



at ~ 30 GeV, similarly to what done in [15], where this model was introduced 
in order to simultaneously fit the CRE fluxes and the positron fraction in a 
two-component scenario (see [TS] for further details). We remark however that 
our qualitative results do not depend on background assumptions. 

If the total flux is given by the contribution of an astrophysical (AP) flux 
and a DM originated one, 4>tot — 4'ap + 4'dm, the degree of anisotropy is given 
by a composition of the AP and the DM intrinsic anisotropics: 



3D V(/)TOT 
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where, again, Sdm = \Sdm\ and Sap 

In a specific scenario, S^ax and Smin are determined by the relative con- 
tribution of (^DM to the total flux. The situation is represented for an energy 
E = 500 GeV by the triangle in Fig. [5]where Sj^ax and Smin si's plotted as a func- 
tion of X = i^DM /4'TOT- The upper side of the blue shaded triangle represents 
Srnax, whilc thc lowcr curve shows Smin- The shaded region between the Smax 
and Smin curves represents the allowed region which the total anisotropy S can 
span. The lower vertex given by the value x — 4>dm / 4>tot = Sap / {Sap + Sdm) 
represents the particular case in which Vc/jdm and Vc/jap are equal and point 
towards opposite directions and, therefore, S = (see Eq. |5|. 

More specifically the left panel of Fig. [s] shows the case of total (AP+DM) 
anisotropy above 500 GeV for a 3 TeV DM, NFW profile, iX^ annihilation 
channel and KRA propagation setup. In this case from Fig. [3] we can see that 
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Figure 5: Allowed region of total CRE dipole anisotropy as a function of the relative DM 
flux in the mixed (astrophysics + DM) scenario. For a given DM flux the allowed range 
for 5 is given by the thick gray segment determined by the intersection of the vertical line, 
corresponding to 4'DM / 4'TOT , with the blue shaded triangle. See text for more details. 



^Hm 1-3 X 10~^ while Sap 5 x 10~^ so that the maximum anisotropy is 
given hy 6 = 5dm achieved when 4>dm ^ (f'AP- This special case is actually 
fairly representative of the general case. From Fig. |3j in fact, it is always true 
that 5dm > 5ap- At the same time, the value 6^'^ ~ 2 x 10~^ at 500 GeV is 
basically model independent as it can be better seen in Fig. |4] 

It is also clear that if the DM and AP fluxes are comparable, (pDM ^ <t>APi 
the maximum anisotropy will be always lower than the maximum anisotropy 
in the case in which DM dominates the flux. Indeed, this intermediate situa- 
tion represents the most likely scenario, since the combined fits of the Pamela 
positron fraction and Fermi CREs already strongly disfavor the case in which 
DM is the dominant flux component, although still allowing a substantial con- 
tribution. 

To check a realistic scenario we consider the case in which DM-I-AP is re- 
quired to not exceed the CRE flux measured by Fermi. In this case we find that 
the CRE flux produced by a 3 TeV DM candidate can have at most a boost 
factor oi B — 470 with respect to the flux produced with the usual thermal 
value of the annihilation cross section {<jav) = 3 x 10~^^ cm'^s"^. In this case 
the DM accounts for a a; = <j)DM/(f'TOT = 82% of the CRE flux at 500 GeV. 
Given the value x = 0.82, the actual allowed range for 5 can be determined by 
the intersection of the vertical line with the triangle as shown in Fig. [5) The 
thick gray segment shows the allowed range for the anisotropy -the actual value 
being eventually determined by the relative directions of '^4>dm and \/<f>AP- A 
detailed fit of the Pamela and Fermi data by various groups [HI H] gives a similar 
or lower DM contribution than the upper limit above, which would give an even 
smaller range of allowed anisotropy 5. The right panel of Fig. [5] show the same 
case but for energy above 100 GeV. Above this energy DM accounts for only 
14 % of the total flux. 

It is clear from the discussion presented above that in any case 5max < 
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max{SDM,SAp), and also that for standard astrophysical background and ac- 
counting for DM substructures, the maximal DM anisotropy is significantly 
larger that the AP anisotropy. Therefore, the maximum DM anisotropy S^'^f 
constitutes an upper limit to the theoretical CRE anisotropy in such a scenario. 
Being S'g'^f very little sensitive to the specific DM scenario and to the details 
of the CRE propagation (see Figs, [s] and [i]), this upper limit is very robust. 

If a positive detection of anisotropy will occur in the future, and the anisotropy 
will be found larger than Sdm, we can then exclude the presence of a substantial 
DM contribution, and therefore we have to demand Sap > 5dm- This would 
point to a scenario dominated by local, discrete astrophysical sources, such as 
pulsars, as the main source of high energy CRE, and would strongly constrain 
the DM contribution to high energy CRE fluxes. However, this argument does 
not exclude that a subdominant contribution from DM annihilation in substruc- 
tures can still be present |4T] (the evaluation of this subdominant contribution 
will depend on the precise characteristics of the model of local astrophysical 
sources and is beyond the scope of the present work). We remark, moreover, 
that in a realistic scenario, as the one discussed above, the maximal theoretical 
anisotropy would be lowered only by about 30% (see Fig. [s]) with respect to the 
theoretical maximum we quote in Fig. |4] 

On the other hand, if only upper limits will be placed and turn out to be 
smaller than the maximum anisotropy, constraints on (j)DM will be possibly 
placed, using e.g. Fig. [s] (although an improved discussion within a three com- 
ponent scenario also including local astrophysical sources will be needed in that 
case.). 

6. Discussion and Conclusions 

Our findings result from a MonteCarlo computation of the local distribution 
of DM substructures and a possible bias of this approach is that we might 
have missed configurations whose probability is less than 1%, in which, e.g., a 
large mass clump emerges isolated and very close to the Earth. This could in 
principle produce an anisotropy larger than what we quote as a "maximum". 
We checked, however, that this configuration cannot produce a high degree of 
anisotropy. Indeed, even in the unlikely case of a 10* Mq clump at 100 pc from 
Earth (whose probability is < 0.1% ^42]) the anisotropy is strongly suppressed 
by the nearly isotropic flux of the much more abundant smaller substructures 
and it is thus always diluted below m&y.{5oM)- This feature makes the DM 
signal intrinsically different from the one expected from pulsars. Indeed, while 
there might be a close-by, isolated pulsar, that can possibly lead to a large 
anisotropy |15j . it is not possible to reproduce this configuration with DM. The 
situation is different also from 7-rays, where this clump would be a quite bright 
point source. 

Another possible caveat is that low mass clumps [nici ^ 10~^ -^0) are so 
abundant that in principle they can be found within 1 pc from Earth, hence 
CREs could reach the Earth before diffusing significantly. Based on their num- 
ber density, we expect to find only a few substructures with mass 10^^ Mq 
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Figure 6: Left: Black solid contours: anisotropy at 500 GeV due to a single clump as function 
of its mass and distance from the Earth when substructures of mass down to 10~® Mq are 
included. Grey dashed contours: expected number of clumps of a given mass closer than 
a given distance to the Earth. Numbers represent the logj^Q of the related quantities. The 
plot is drawn for 3 TeV DM fully annihilating into muon pairs, assuming a NFW profile and 
KRA propagation setup. Right: same as left panel, but when no substructures other than the 
clump itself (and no main Halo) are included. The white dashed curves represent the logio of 
the boost factor (with respect to thermal cross section) required for the CRE fluxes to be in 
agreement with the observed fluxes at 500 GeV. 

within 1 pc from Earth. These clumps would look more like point-like sources 
of e+e~ rather than like a dipole. Even in this case, however, their point-like 
flux both in e'^e~ and 7-rays would be several orders of magnitude below the 
Fermi sensitivity. 

These points are illustrated in the left panel of Fig. [6] where the contours 
show the anisotropy at 500 GeV of a single clump as function of its mass and 
distance from the Earth. The plot represents the same case assumed in drawing 
Fig. [5} 3 TeV DM fully annihilating into muon pairs and assuming NFW profile 
and KRA propagation setup. The anisotropy is calculated as in Eq. [4] with at 
the numerator the gradient of the flux of the clump and at the denominator 
the total flux from all the clumps and the smooth Halo so to represent the 
effective contribution of the clump to the total anisotropy. As it can be seen, 
the anisotropy has a plateau at about 10~^ (corresponding to the case where 
the observer is well inside a very massive clump) which does not exceed the 
maximum possible anisotropy at the same energy. 

A possible exception to the above scenario is the extreme case in which 
a massive DM clump is the "only" relevant CRE source. In fact, any other 
contribution from other clumps or the smooth halo would dilute the anisotropy 
of this single clump. We note here that in this scenario, in order for the halo 
contribution to be negligible with respect to the one of the single clump, one 
needs to invoke a strong suppression of the annihilation cross section of the DM 
in the halo much below the thermal value. Moreover, clumps of at least 10^ Mq 
are known experimentally to exist in the form of dwarf galaxies |43| and they 
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would share the same boost factor as the hypothetical nearby clump, so that 
their presence would dilute the anisotropy of the single clump anyway. This 
case is discussed here only because of its extreme geometry and is analyzed in 
the right panel of Fig. [6j The anisotropy is again calculated as in Eq. |4] but 
this time with at the numerator the gradient of the flux of the clump and at 
the denominator the flux of the clump itself. The needed boost factor with 
respect to the thermal cross section in order for the clump to produce all the 
observed CRE flux at 500 GeV is shown as the white, dotted curves. This 
gives the maximal possible anisotropy and, as we expected, it can exceed the 
upper limits of Fig. |4j being possibly as large as 10~^. However, the only 
clump configurations which would give such a high anisotropy lie in a region of 
parameter space (small mass or large distance) where unrealistically high boost 
factors (of the order of 10^ or higher) are required. If we restrict the allowed 
region to the more reasonable case of a nearby massive clump (lower right corner 
of the plot), again the maximal anisotropy does not exceed ~ 10~^ since the 
anisotropy is somewhat reduced by the fact that the observer is well inside the 
clump and sees the structure of the clump. We notice that large boost factors 
(e.g. via Sommerfeld enhancement) are strongly constrained by analyses of the 
CMB distortion during the recombination epoch [44j|45l|46]. Further constraints 
on the boost factor also come for DM annihilation in the core of the Earth [57] . 
Also, in the case of a small clump, besides unrealistically high boost factors being 
required, it would be not justified anymore to not include the contribution of 
the other equally or more massive clumps, falling back to a configuration of low 
anisotropy similar to the one of Fig. |6] (left). For both cases it should be further 
stressed that, although above ~ 200 GeV the positron fraction and thus the 
possible DM fraction is unconstrained, it is unlikely that DM constitutes 100% 
of the CRE flux. Considering the contribution of the astrophysical component 
will lower the overall anisotropy of the clump although the precise decrease will 
be dependent on the particular AP model employed. 

Another remark concerns the density profiles we considered. While iV-body 
simulations suggest spiked halo and subhalo matter density profiles, astrophysi- 
cal observations of many dwarf spiral galaxies point to a shallower, Burkert-like 
density profile [15]. Our results are stable under the relevant change from a 
spiked to a cored profile. Indeed, high energy CREs arriving at Earth do not 
carry information on the DM distribution in the galactic center, as they propa- 
gate only a few kpc in the interstellar medium. The anisotropy is not sensitive 
to the internal concentration of the subhaloes as well, because diffusion over kpc 
scales smooths out the effect of a possible cusped over-density region. For the 
same reason, in the case of decaying DM we find similar results as in the case 
of annihilating DM. 

We also neglected the effects of a possible proper motion of substructures. 
Indeed, as it was pointed out in [IS] for the case of an isolated substructure, a 
dynamical treatment would lead to a slightly enhanced dipole anisotropy only 
for sources moving towards the Solar System. However, while this effect can be 
relevant for a single clump, it is expected to average away for a population of 
clumps as considered here. 
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Finally, a word of caution must be said about our choice of diffusion models. 
It might in fact be that the local anisotropy observed on Earth can be affected 
by local magnetic turbulence, which would break the assumption of isotropic 
and spatially uniform diffusion we used in this work. Hints in this direction 
may come from the observations of dipole anisotropies and on anisotropics on 
angular scales of the order of 10° 30° in the CR hadronic component at 
energies > 10 TeV (SOI ED E3 EH ES], as discussed in several works [53 
E3 EH ESI EQI m] . Anisotropies typically increase as a function of energy on 
account of the increasing gyro-radius, and thus the intensity of the anisotropies 
observed above ^ 10 TeV will be correspondingly decreased when rescaled to 
our energy range (^100 GeV - 1 TeV). On the other hand leptons in this energy 
range have a much smaller horizon with respect to hadrons (due to their shorter 
propagation length) and this typically can increase the anisotropies, since local 
effects are more important. Thus, it would be difficult to understand which is 
the dominant effect and to assess precisely the effects of local turbulence on our 
results. We remark, however, that the observed hadronic dipolar anisotropy 
seems generally in agreement with the hypothesis of isotropic and homogeneous 
diffusion [TB], while local magnetic turbulence seems to be required mainly to 
explain the anisotropy at higher multipoles [611 . Since we are considering in this 
work only dipolar anisotropies, our results are likely less affected by the above 
effects. 

In summary, barring the above caveats, we demonstrate that our results on 
DM anisotropy are robust with respect to several choices of propagation setup 
and of DM spatial distribution and particle model, and we thus propose to use 
them as a criterion to reject or at least disfavor a DM dominated scenario in 
the case of detection of a large anisotropy in high energy CREs. 
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Appendix A. Solution to the diflfusion equation 

Assuming spatially uniform D and b and looking for a stationary solution, 
Eq. ([3| reduces to the form: 



D{E) A 



dn 
dE 



d_ 
dE 



b{E) 



dn 
dE 



Q{E,^) . 



(A.l) 
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We look for the Green function G(x, E, x', E') so that the solution of the equa- 
tion can be written as 



dn 
dE 



iE,x)^J dE' J d^x' G{x,E,x',E')Q{E',x') . (A.2) 



Following im IMl US] , the Green function Gp for the free case (i.e. with no 
boundary conditions) can be written as 

Gpi.,E,.',E')^J-(-^X^\^j>(~\^] , (A.3) 



b{E) yAnAr J ^ V 4A 
where we define 

r^' Die) 

Ax = X - x' and Ar = / vT^de . (A.4) 



6(6) 



Therefore, 



where Pe// is such that 



Peff 



In the case of a clump 

p(x') = Po/(|x'-x,,|/ro) , (A.7) 

with Xci being the position of the center of the clump. 

By changing variables to spherical coordinates local to the clump's center, 
and using ^ = |x' — Xc;|/ro, we end up with[^ 

(A.8) 

An interesting limit is given by the situation in which the observer is well 
outside the clump. In this approximation we have |Axp = |x— x'p ~ |(x— Xci)^ 
and the above expression simplifies as 



^For a NFW profile x_2 = 1, while for a Burkert profile x_2 = (1 - ^26/27)1/^ + (1 + 
^26/27)1/^ ~ 1.52. 
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Table B.l: Density parameters (unbiased case) 





rah 


POh 


POcl 


Ch 


A 




kpc 


GeV c^^cm^^ 


GcV c^^cm^^ 




Mq^ kpC-3 


NFW 


21.7 


0.132 


0.151 


10.9 


1.08 X 10^ 


Burkert 


13.5 


0.404 


0.462 


11.5 


3.30 X 10^ 



where s = |x — Xc/|. 

The decaying DM case is completely analogous. We simply have to substitute 
((TAw)/2 with r and p^^^ with a peff defined in terms of /3(x'). 

Appendix B. Distribution of substructures 

We assume that the number density of subhaloes scales with mass and po- 
sition according to 

^'-(to,z,x)=A gcii-] , (B.l) 



dnici \Mq J \ro 

where gdir/ro) is a dimensionless function parameterizing the spatial distribu- 
tion of the substructures; tq is a scale parameter for the entire DM distribution; 
A is a dimensional normalization constant. The effect of tidal disruption of 
clumps near the galactic center is included a posteriori in our Monte Carlo 
according to the Roche criterion (see e.g. [33]), but not in the following cal- 
culations which are normalized to clumps with mass < IO^Mq and are thus 
unaffected by the small number of tidally disrupted clumps. Similarly we only 
use a = 2, because with a minimum mass scale of 1O~^M0, a mass index of 2.0 
or 1.9 produces only a minor change in the results. 



From Eq. (B.l ), the total mass and number of DM clumps of mass between 



mi and TO2 can be written as 



/I*m2 d/i 
d^x / m,ci-—^{mci,yi)Amci (B.2) 



dn 



N{mi,m2) = / d-^x / ^p^(mc;, x)dmc; . (B.3) 

According to Via Lactea II [30] we impose the condition that the mass con- 
tributed by clumps of mass between lO^M© and lO-^^M© is 10% of the total 
mass of the Milky Way, M@. From the above formulae we can then see that: 

16 

Ma - M(1O-^M0, IQIOMq) ^ — M@ ~ 53.3 %M@ 

oU 

1 0^ A/f 

7V,z = 7V(10-6Mo,101°Mq) = — — -?^^2.90xl0i^ (B.4) 

30 m 10 
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Table B.2: Density parameters (antibiased case) 





roh 


Po 


n 


A 




kpc 


GeV c^^cm^'^ 


kpc 


Mq ^kpC-3 


NFW 


21.7 


0.284 


61.5 


2.03 X 10^ 


Burkert 


13.5 


0.866 


54.4 


6.19 X 10^ 



The DM density is the sum of its two components, the smooth halo (h) and 
he clumpy one (cl): 

Ptot{r) = Ph{r) + Pci{r)- 

The spatial distribution is less known, and two different hypotheses have been 
proposed till now, an unbiased distribution which assumes that the subhaloes 
distribution follows the same radial profile of the main halo, and an antibiased 
distribution which assumes that the two are anti-correlated. We will consider 
both cases separately. 

In the case of unbiased distribution for the substructures we have that ph 
and Pel have the same spatial dependence: Pi{r) = Poifir/ro), in particular 
they share the same scale parameter tq. The three parameters tq, poh and 
Poci describing the DM distribution can be determined using the three known 
quantities M@, the total mass of the Milky way, ps, the local value of the DM 



density, and the total mass in substructures derived in Eq. (B.4). Table B.l 
shows the results obtained assuming M@ = 1.49 x 10^^ Af©, Rs = 8.28 kpc, 
PS = 0.389 GeV c"^ ^66j and two different choices for the DM density. For 
M@, Rs and ps we use the values obtained from [66] assuming a NFW density 
profile, because they are determined with a slightly better accuracy. The values 
obtained under an Einasto hypothesis are compatible with these within one 
standard deviation 

For the anti-biased case we follow the approach of [33] in which the original 
antibiased (cored) distribution of subhaloes, initially proposed in [67l [68], is 
modified in order to be consistent with a given overall density profile of the 
Galaxy. The two contributions PhiT) and Pdif) to ptot{f) are written in the 
form ^ 

Ph{r) = — -i—ptotir) 

l + r/n 

Pci(r) = — ]—ptot(r) 

l + r/rt 

in terms of the bias radius ri,, determined by inverting the equation 

Ph{r)4TTr^dr = = M@ - 



See [33] for more details, ro and pa are the parameters defining ptot = Pof{r/ro). 
The results are shown in table IB. 21 

The last piece of information required to derive the annihilation signal from 
the clumps is the DM distribution within the clumps themselves. We will assume 
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Figure C.7: Boost factor given by dark matter substructures to the electron and positron 
fluxes. The sohd line shows the case of unbiased DM substructure distribution considered in 
the paper. The shaded region represents the fluctuations from different realizations of the 
substructures. The dashed black line corresponds to the antibiased distribution used also in 
| 33| which is instead shown with the red dashed line. 



that the clumps foUow the same mass profile as the main halo, but with their 
own parameters r^ ci and po,c/ replacing the ones associated to the halo. We 
parametrize the internal concentration of the subhaloes like in |33i ( Via Lactea 
II case). In particular, since we deal with the smallest clumps, we fit it to a 
simple power law in the mass range 10^^ 10^ M©: 

Ccl = Co (md/Mg)"^ 

with Co = 102.8 and /? = 0.0331. 



Appendix C. Mean Electron and positron flux and gradient 

The mean flux from all the substructures can be expressed as: 

due drici , „. 



b{E) An Je 2 m\ dE' 

A similar expression holds for the mean gradient. The rather complicate equa- 



tion above actually quite simplifies in the "far clump" approximation (see Ap- 



pendix A) where the integration over the clump positions and masses can be 



performed separately and evaluated in a straightforward way. 
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We show in Fig. |C.7| the boost factor of DM substructures to the electron 
and positron fluxes and we compare it with the same quantity computed in 
[33] for a similar DM distribution. The solid line shows the case of unbiased 
DM substructure distribution considered in the paper and the shaded region 
represent the fluctuations from different realizations of the substructures. In 
order to compare our results with those shown in [33] we also show the case 
of antibiased substructures' distribution. As it is clear from Fig. |C.7[ passing 
from an unbiased to an antibiased model reduces the relative contribution by 
roughly a factor of 3, in very good agreement with the results of [33]. The 
residual difference is due to some differences in the global parameters we adopt 
to describe the dark matter distribution. These differences, however, are not 
relevant when computing the anisotropy. 
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